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Abstract 

By means of the perturbative renormalization group method, we 
study a long-time behaviour of some symplectic discrete maps near 
elliptic and hyperbolic fixed points. It is shown that a naive renor- 
malization group (RG) map breaks the symplectic symmetry and fails 
to describe a long-time behaviour. In order to preserve the symplec- 
tic symmetry, we present a regularization procedure, which gives a 
regularized symplectic RG map describing an approximate long-time 
behaviour succesfully. 



1 Introduction 

There has been a long history to study an asymptotic solution of Hamil- 
tonian flows by means of singular perturbation methods such as the averaging 
method and the method of multiple time-scales. A Hamiltonian flow can be 
reduced to a symplectic discrete map called the Poincare map, which has the 
lower dimension than the original flow and is, therefore, extensively studied. 
However, neither the averaging method nor the method of multiple time- 
scales may be immediately applicable to symplectic maps. The perturbative 
renormalization group (RG) method developed recently may be a useful tool 
to tackle asymptotic behaviours of discrete maps as well as flows. The original 
RG method is an asymptotic singular perturbation technique developed for 
diffential equations Secular or divergent terms of perturbation solutions 
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of differential equations are removed by renormalizing integral constants of 
the lowest order solution. The RG method is reformulated on the basis of a 
naive renormalization transformation and the Lie group 0. This reformu- 
lated RG method based on the Lie group is easy to apply to discrete systems, 
by which asymptotic expansions of unstable manifolds of some chaotic dis- 
crete systems are obtained 0. The extension of the RG method to discrete 
symplectic systems is not trivial because the symplectic struture is easily 
broken in naive renormalization group equations (maps) as shown in this pa- 
per, while the application of the RG method to Hamiltonian flows does not 
cause such a problem as the broken symplectic symmetry except a special 
case 0. The application of the RG method to some non-symplectic discrete 
systems was tried in the framework of the envelope method Howevever, 
the method, if applied to a symplectic map, would give only a naive RG map 
which breaks the symplectic symmetry. 

The main purpose of the present paper is to present an RG procedure to 
preserve the symplectic structure in RG maps and to obtain symplectic RG 
maps. In this paper, this procedure is called a regularized RG procedure, 
which consists of the following two steps. First, using the reformulated RG 
method , we get a naive RG map near elhptic and hyperbolic fixed points 
of some symplectic discrete systems. The naive map preserves the symplectic 
symmetry only approximately and fails to describe a long-time behaviour of 
the original maps. Second, in order to recover the symplectic symmetry, we 
introduce a process of "exponentiation" which yields a symplectic RG map. 
This process and a symplectic RG map thus obtained are called,respectively, 
regularization of an RG map and a regularized RG map. 
In section 2, a long-time behaviour of a simple linear map is analyzed in 
order to elucidate the broken symplectic symmetry in a naive RG map and 
our regularization process. In section 3 , regularized (symplectic) RG maps 
are obtained near elliptic and hyperbolic fixed points of a two-dimensional 
nonlinear symplectic map. In section 4, the regularized RG procedure is 
applied to a four-dimensional symplectic map, of which elliptic-elliptic fixed 
point has an irrational frequency ratio. 
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2 Linear symplectic map 

It may be instructive to analyze a linear symplectic map, which is exactly 
solvable, 3 {xn,yn) ^ {Xn+l,yn+l) ■ 

Xn+1 Xn ~\~ Un+li 

Un+l = Vn- Jx„ + 2eJXn, (l) 

where n E Z and e is a small parameter. The map (|I]) has a elliptic fixed 
point at the origin (0,0) for < J < 2. Eliminating y from (|I|) ,we obtain a 
second order difference equation: 

LXn = Xn+l - 2 C0S(6')X„ + Xn-l = e2JXn, (2) 

where cos 6' = 1 — J/2 and < J < 2 is assumed. The linear map (§) has 
the following exact solution x^ : 



x^ 



Aexp arccos(cos6' + eJ)n^ + c.c. 



Aexp 



sm 6 2 sm 6 V sm 6 J 



+ C.C., (3) 



where A(g G) is a complex "integral" constant and c.c. stands for the com- 
plex conjugate to the preceding terms. 

Let us derive an asymptotic solution of the map (|I|) for small e by means of 
the RG method. Substituting the expansion : 

= + 6X« + e^xl^) + • • • , (4) 

into Eq.(|^), we have 

= 0, Lx^y = 2 Jxi°) , Lxl^) = 2 Jxi^) ■ • • . 



and 



x^^ = A exp{i6n) + c.c, 

x^}-^ = — — -n exp(z6'?T,) + c.c. 
i sin 9 

-J'^A / o .cos 6* 



(2\ -J A ( 2 .cosy \ , 

"Ai = ^ . 9 .1 — -n] exphtin) + c.c, 
2sm^/9V sm6' / 
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where C) is a complex constant. To remove secular terms (oc n, n^), we 
introduce a renormalization transformation A hh^ A{n) [Q: 



JA 



A{n) = A + e . n + e 



-J^A 



.cos 6 



tsmO 2 sm 6\ sm6 ' 



A discrete version of the RG equation is just the first order difference equation 
of A{n), whose local solution is given by Eq.(H). From Eq.(|), we have 



A{n + 1)-A{n) 



— e 



sin d 2 sin^ Q \ 



f2n + l + z^^ )A + 0(e3), (6) 



sin Q 



where A should be expressed in terms of Ain). This is done by taking the 
inversion of the renormalization transformation (0) iteratively. 



A 



l + ^e^ + 0(e2))A(n). 
sm y 



(7) 



Substituting (|^) into we obtain the following RG equation (RG map) up 
to 0{^). 



A{n\\) 



1 + 



—ieJ 1 f—ieJ\'^ 
+ 



sin 9 2! V sin 9 



— le 



; COS 9 

2 sin^ 9 



A{n) + 0(6=^ 



of which solution is 

I(n) = 



—ieJ 1 f—ieJ\'^ 



sin 9 2! \ sin 9 



le 



, cos 9 
2 sin^ 9 



+ 0{e')] A{Q). (9) 



On the other hand, from Eq. (j^), we have A{n) exactly as 



A{n) = A{0) exp 



-J 



cos9 / J \2 



sin^ 



2 sin 9 V sin 9 



+ 



n 



(10) 



Notice that |y4p is an exact constant of motion while it is merely an approx- 
imate conserved quantity of the (trancated) RG map (||). The symplectic 
structure is also not exactly preserved in the RG map, that is, for the trun- 
cated RG map (||) up to 0{e''), we have 



dA{n + 1) A dA*{n + 1) - dA{n) A dA*{n) = 0{e 



(fc+l)^ 



7^0, 
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where = 1,2, A* is complex conjugate to A and should also be a 
canonical conjugate to A ^j. Although this fault of the RG map vanishes in 
the limit of small e, it is intolerable as shown in Fig.l where we depict a long 
time behaviour of the solution for small but finite e. In fact, the truncated 
RG map becomes a system with small but finite dissipation . In order to 
remedy a fault like this, we take advantage of a crucial observation that the 
coefficient of A{n) in (H) can be modified as 

—ieJ 1 /— ieJx^ 2^^cos^^ 



sine 2\\sm9J 2sm^e 

exp ^ - le^ ^ + O . 

^Vsin^ 2sin^e/ ^ ^ 



Using this modified coefficient, the symplectic symmetry in the truncated 
RG map is recovered as 

/l(n + l)=exp TT-^e^ (H) 

which also has an exact conserved quntity |y4|. This process is nothing but 
"exponentiation" of the coefficient of A{n) and is called regularization of an 



RG map. It is easy to see that the solution of the regularized RG map (|TT 



coincides with the asymptotic expansion of the exact solution (|10D and de- 
scribes a long time behaviour up to n ~ 0{e~^) in the present approximation. 
In Fig.(l), trajectories constructed from the naive RG map and the reg- 
ularized RG map (|TT]) are depicted to be compared to an "exact" trajectry 
of the original map ([l|) 
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X 



Figure 1: Trajectries constructed from the naive RG map (broken line), the 
regularized RG map(solid line) and an "exact" trajectry of the original map 
([l|) (crosses) for J = 0.2, e = 0.25 and an initial condition {xq = 1.0, i/q = 0.0). 



3 Two-dimensional Non-linear Symplectic Map 

3.1 Elliptic Fixed Point 

Let us analyze a weakly non-linear symplectic map 

3 {Xn,yn) ^ {Xn+l,yn+l) ■ 

Xn+1 -\- Xjin+l^ 

Vn+l = Vn- JXn + '^^eJxl, (12) 

or 

Lxrr = e2Jxl, (13) 

where e is a small parameter, L is defined in (^ and < J < 2 is assumed 
so that Eq. ([T3|) has a elliptic fixed point at the origin (0,0). Expanding x„ 
as a power series of e 

x„ = xW+ea;« + e2xi2) + ..., (14) 

we have 

La;W=0, Lx«=2Jx^, Lxf = aJx^x^, ■ • ■ , 
and solutions of the perturbed equations to 0{e^) are given by 
x^^-* = Aexp{i6n) + C.C., 

x^^^ = ai exp{3i6n) + iaiR\A\'^Anexp{i6n) + C.C., 

{ COS Q \ 

x)^ = iaiaiii\A\^Anexp(i6n) ^lAI'^Al + i n ) exp(i6n) 

2 \ sin6' / 

\ A\i sin(3^) \ 

-l-Szaiaiij A A \ n — i expfSit^n) 

V cos W — cos 6 



+Qai\A?' A^ ex-pi'iien) H "^^^^^ — -exp(526') + c.c, 

cos bO — cos 6 



where 



_ J _ -3J 

C^l = 7^ OllR = 



COS 3^^ — COS 9 sin 9 

and A{& C) is a complex integral constant. In order to remove secular terms 
in the coefficient of the fundamental harmonic (exp(i^n)), we introduce a 
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renormalization transformation A i-^ A{n): 
A{n) = A + e{iaiR\A\'^An) 

+,^( zdli\A\'A(n' + z^n) + za,a,R\A\'An] . (15) 
\ 2 V sinb* / / 

Following the same procedure as that in the preceding section, we derive a 
naive RG map from Eq. ([l5|) : 

A{n + 1) - A{n) = eiaiR\A{n)\'^A{n) 

— a?R,T/ M/i T/ X .cos^^ 



+iaiam\A{n)\^A{n)), (16) 



which breaks the symplectic symmetry and does not have a constant of mo- 
tion. To recover the symplectic symmetry, we regularize the naive RG map 
([l6| ) by "exponentiating" the coefficient of A{n), that is. 



A{n + 1) = A{n)exp yie\A{n)faiR 

+^e2|I(n)|^(^^a?H + "laii?)) , (17) 

which of course agrees with Eg. ([TB|) to O(e^). It is easy to see that Eq.(p!7D 
has the symplectic symmetry 

dA{n + 1) A dA*{n + 1) = dA{n) A dA*{n), 

and a constant of motion |y4|. Introducing the polar expression A{n) = 
\A{n) \ exp(i0(n)), the regularized RG map is reduced to a simple phase equa- 
tion given by 

0(n + l) = 0(n)+ (^e|I(0)|V/? 

+e^l^(0)r(^^a?^ + aiaiK)), (18) 



Thus, the regularized RG map (jl^) is exactly solvable in contrast to the 
naive RG map (0) which do not have a conserved quantity. It is noted that 
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secular coefficients of the third harmonic (exp(3i^n)) are also summed up to 
give a renormalized coefficient A^ln) as: 

+e' (3za^a,n\A\'A' (n - if J + 6al\A\'AA .(19) 
\ V cos St' — cos t^/ / 

Substituting an iterative expression of A in terms of A{n) obtained from 
Eq.(|l^) into Eg. (]T9|) , we can eliminate secular terms in the same way as the 
case of differential equations 0. 

As an application of the result, we obtain an approximate but analytical 
expression of the rotaion number near the elliptic fixed point. The definition 
of the rotation number is 



N 



n=0 



where 0„ is an angle between the vectors {xn,yn) and {xn+i,yn+i) 0- Ne- 
glecting effects of the higher hamonics, an approximate expression of p{xq, yo) 
of the present system (KT) is gven as following : 



2TTpixo,yo) ^ e + (^e\A{0)\'aiR + e^\A{0)\'(^^^alR + aiaiR^^ +0(0), 

(20) 

where initial values {xo,yo) are related to |v4(0)| and 0(0). In Fig. (3), the 
rotaion number p(xo,?/o) given by (^) is depicted for initial values on the 
half line < xq < 0.8, yo = shown in Fig. (2). The result agrees well 
with the "exact" rotation number obtained by numerical calculations of the 
original map (|1^) until a new resonance appears near xq ~ 0.6. 
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X 

Figure 2: A phase portrait of the map (|12]) for J = 1.2, e = 0.25. The sohd 
half line represents initial phase points < xq < 0.8, yo = of the rotation 
number p{xo, yo). 
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Figure 3: The rotation number p{xo,yo) versus initial phase points (xo,0) 
for J = 1.2, e = 0.25. The sohd curve represents the rotation number 
given by the regularized RG map (PP]), while the dashed curve is obtained 
by numerical calculations of the original map (p^. 
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3.2 Hyperbolic Fixed Point 



In the case J < in ([T^), the origin (0,0) is a hyperbohc fixed point of 
T^ ) and Eq. ([T3D becomes 

Lxn = e2Jxl, (21) 



where 



Lxn = Xn+i — 2 cosh.{9)xn + Xn-i, cosh6' = l — J/2 (J < 0), 

and K = exp(^) is one of the eigenvalues of the hnearized equation at the 
origin, that is 



^(2-J + V-4J + J2 j , or K = ^(2-J- V-4 J + J2 
Expanding the (|14|), we get the following perturbed solution 



K 

2 



T^^> = 2 71— t \ = \ - riK A ~ nK 

13 /l2 , , 43 a2 



^(2) 



+ 12 J2 ( ^ ^"^f f + a{K)n) + ^^^f—-^ i^-" + n.r., 

\2{K-^-Ky\ ^ ^ ) {K'^-K)D{K)) 



where 



D{K) = K-^ - {K + K-^) + K-^ = D{K-^) e R, 
a{K) ^ j^^^ = -a{K-') e R. 

and A- are real "integral" constants and n.r. denotes non-resonant terms 
which are proportional to K^"', K~^"', K^"', K~^^. Introducing a renormaliza- 
tion transformation A^ h-* ^+(^) and A_ ^ A^(n) to remove secular terms 
in the coefficient of fC", we have a set of naive RG equations: 

T . ^ T / ^ 3li(n)I_(r2) 
AJn + l) = AJn)-2eJ- +^ ' ^ ' 



K-K-^ 

+12eV^ 



3At(n)Ai(n) 
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A_{n + 1) = A4n)-2eJ 



{K - K-^)D{K))' 

3Al{n)A+{n) 



(22) 



K-^-K 



f 3Al{n)Al{n) f^ , 

(,2(^^r3^(;i + «W 



_4i(ri)A^(n)_\ 
+ ^^^^ 

The set of Eqs. (^21) and (^) do not have the symplectic symmetry but is 
regularized by the folowing exponentiation procedure. 

7 / T / N f 6jI+(n)I_(n) 

A+{n + l) = A+{n)exp{e +^ ^ ^ ^ 



V ^ ^ 2(ir-i - K)^ {K-^ - K)D{K) ) ) ' 

This regularized RG map has a constant of motion 

P = = I+(0)I_(0), 

and is also exactly solvable. In terms of P, a general solution of Eqs. 
and (p5D is given as 



A4n) = A+{0)exp(Q{P;e)n 
A^{n) = l_(0)exp (^-Q(P;e)n 



where Q{P; e) is a polynomial of P: 



\2{K - K-^Y {K - K-^)D{K) 
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4 Four- dimensional Symplectic Map 



Let us consider a coupled map of two symplectic maps (|T2|) 

dF dF 



X. 



n+1 



dyn+1 ' 
OF 



n+1 



dXn' 

OF 
dx' '' 



(26) 
(27) 



n+1 + Vn+lj 



where F is a gerenating function 

^ [^ni Un+li ^ni Hn+l) 

+ + J'^^) - 2e(j^ + J'^^ + axlx':). 

The origin (0, 0, 0, 0) is an elliptic fixed point of (^ and (g^) for < J < 2 
and < J < 2. Eliminating the y and y from Eqs. (p6D and (|27|) , we obtain 
coupled second order difference equations: 

Xn+i - 2 cos(6')x„ + a;„_i = 2e( Jx^ + 2aXnx'^), (28) 

x^+i - 2 cos(^')a;L + x'^-i = MJ'^n + "^axlx'j, (29) 

where cos ^ = 1 - J/2 , cos 0' = 1 - J /2. Now, we concentrate on the case 
that the ratio of frequencies of the leading order solution of ( pHj) and (|29D 
,i.e. 6/6' is an irrational number. Then, straightforward calculations yield 
perturbation solutions of (p8|) and ( pQ]) as 

= Aexp(i6'n) + c.c, 



XI 



X, 



'(0) 



i? exp(i6' n) + c.c. 



a; 



(1) 



X 



'(1) 



— z 

sin ^ 
—i 

sin 



(3J|ApA + 4:a\B\'^A)nexp{in6) + c.c. + n.r.. 



(30) 



j{3j'\B\^B + 4:a\A\^B)nexp{in6') + c.c. + n.r., (31) 



where A G C and B E C are integral constants. In order to remove secular 
terms in (|30|) and (31), we introduce a RG transformation A i— > y4(n),i? i— »• 
B{n): 

A{n) = A + e^iSJlAl"^ + 4a\B\^)An, 



Bin) = B + e 



sin 6 
—i 

sin 6 



j{3j'\B\^ + 4:a\A\^)Bn, 
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from which a naive RG map is obtained as 

A{n + 1) = A{n)\l + e^{3J\A\^ + Aa\B\^)], 
I sm u J 

B{n + 1) = B{n)\l + e^^{3j'\B\^ + 4a\A\^)]. 

I sm t/ J 

To regularize this map, we scale the renormalized variables as 

_ A{n) _ B{n) 

V sm V sm 6 

where sin^ > and sin^' > 0. Finally we "exponentiate" the resultant naive 
RG map as 

a{n + l) = a{n)exp\e{-t)(3J^^\ain)\^ + 4:a\P{n)\^)], (32) 

I V sm ' ) 

/?(n + l) = /?(n)exp{e(-0(3j'^|/5(n)|2 + 4a|a(n)p)}, (33) 

It is easy to see Eqs.(^) and ( |55| ) have the sympletic symmetry 
(ia(n+l)A(ia*(n+l)+(i/5(n+l)A(i/3*(n+l) = da{n)Ada*{n)+d/3{n)AdP*{n), 
and two constants of motion |a;p and 

5 conclusion 

We present a regularized RG procedure to preserve the symplectic struc- 
ture in RG maps near elliptic and hyperbolic fixed points of some symplectic 
discrete systems. The regularization is accomplished by simple exponentia- 
tion of a naive RG map and gives a symplectic RG map, which succesfuUy 
describes a long-time asymptotic behaviour of the original system. Without 
such regularization, a naive (truncated) RG map fails to describe a long-time 
behaviour especially near elliptic fixed points. Furthermore, all regularized 
RG maps obtained here have constants of motion and can be exactly solvable, 
while non-regularized (nonlinear) RG maps have no conserved quantities and 
are not solvable analytically. As an application of the regularized RG map, we 
give an approximate but analytical exprssion of the rotation number around 
an elliptic fixed point. 
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It is easy to see that the present regularization process is also apphcable to 
general two-dimensional symplectic maps and higher dimensonal symplectic 
maps, which have an elliptic fixed point with incommensurate frequencies. 
However, construction of regularized RG maps of general higher dimensional 
symplectic maps is still an open problem to be studied in future. 
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